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Selftrapping has been traditionally studied on the assumption that quasipar- 
^ , tides interact with harmonic phonons and that this interaction is linear in 

'^ ' the displacement of the phonon. To complement recent semiclassical studies 

o 
o 



of anharmonicity and nonlinearity in this context, we present below a fully 
quantum mechanical analysis of a two-site system, where the oscillator is de- 



H ' scribed by a tunably anharmonic potential, with a square well with infinite 

walls and the harmonic potential as its extreme limits, and wherein the inter- 
action is nonlinear in the oscillator displacement. We find that even highly 
anharmonic polarons behave similar to their harmonic counterparts in that 
selftrapping is preserved for long times in the limit of strong coupling, and 
that the polaronic tunneling time scale depends exponentially on the polaron 
binding energy. Further, in agreement, with earlier results related to har- 
monic polarons, the semiclassical approximation agrees with the full quantum 
result in the massive oscillator limit of small oscillator frequency and strong 



quasiparticle-oscillator coupling. 
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I. INTRODUCTION 

Recent work by Grigolini and coUaboratorsa, and by Salkola and the present authorsafl, 
have uncovered subtle features associated with selftrapping of quasiparticles in interaction 
with vibrations. The vibrations considered in all those analyses have been harmonic. The 
question of how polaron dynamics and selftrapping are affected by anharmonicities in the 
vibrations was raised by KenkreP several years ago at the level of the discrete nonlinear 
Schrodinger equation (DNLSE) and analyzed by Kenkre and collaborators in the context 
of rotational polaronso'Ei, exponential saturation&Q, and general considerations^. Since the 
validity of the DNLSE has been called into question by recent consider ationsETQ, it is impor- 
tant to examine the issue of what polarons, or selftrapping, owe to harmonic features from a 
starting point which is fully quantum. The present paper is devoted to such an examination 
for a two-site system. 

We will focus here on confined systems rather than on periodic systems such as those 
which may lead to rotational polaronsQiH. In a certain sense, the most anharmonic potential 
conceivable is that which corresponds to a box with infinitely high walls as it corresponds to 
a harmonic piece with vanishing frequency throughout the interior of the box but one with 
infinite frequency at the wall. We choose the symmetric Poschl- Teller potential given by 

Vpt{x) = UotSiia^{ax), (1) 

because it allows continuous transition between the harmonic oscillator and the box limits 
and because it can be treated analytically with ease. In Eq. (P f/o, a are constants defining, 
respectively, the strength and confining region of the potential. The potential becomes 
infinitely steep at x = ±7r/2a. By rewriting the strength of the potential Uq as A(A — 
l)fi^a? /2m where m is the mass of the particle, we introduce the parameter A which describes 
the departure of the potential between the box and harmonic oscillator limits. In the limit 
A ^ 1, the Poschl- Teller potential becomes the infinite square well of width vr/a. In the 
opposite limit A — > oo, a -^ 0, \o? remaining constant (and finite), one recovers the harmonic 
oscillator potential. The eigenenergies of the Poschl- Teller potential (P are given bya 



En = ^^{n^ + 2nX + X), n = 0,1,2,... (2) 

2m 

and the corresponding eigenf unctions are 

(j)n{x) = (x|0„) = Nncos^^"^ axP^^^^^y^i^inax), (3) 

where -P^ (t) are the associated Legendre functions, with A^„ = ( " pfa+y — ) 

II. ANHARMONIC POLARON - STATIONARY ASPECTS 

Consider a two-site system consisting of a quasiparticle, hke an electron or an exciton, 
whose inter-site hopping is described by a matrix element of strength V. The quasiparticle 
also strongly interacts with a vibrational mode between the two sites. This vibrational mode 
is described by the Poschl- Teller potential (|l]). In the harmonic case the usual interaction is 
linear in the vibrational amplitude and consequently connects nearest neighbour eigenstates 
of the oscillator. These two features are distinct from each other. In developing a scheme 
for analyzing effects of a genealization to anharmonic situations, we must maintain either 
one or the other of the two features. We have studied both cases. Here we present results 
of maintaining the second feature, viz., an interaction which joins nearest neighbour energy 
eigenstates. As Nieto and Simmons^ and Crawford and Vrcsaytl have pointed out in a dif- 
ferent context, a sinusoidal interaction posseses this feature for the Poschl- Teller potential, 
and also reduces to the linear form in the harmonic limit. We thus take the full Hamiltonian 
of our system to be 

2(A+l/2) 
where 

.o=^(A + i) (6) 

is the difference between the energies of the first excited state and the ground state of 
the Poschl- Teller potential, z is the dimensionless oscillator coordinate ax, and g is the 
quasiparticle-oscillator coupling constant. Here and henceforth we put h = 1 for simplicity. 



(vf/ + A(A - 1) tan^ z) + ujogjfx + -jp smz + Vf, (4) 



The operators p, f are the operators describing the quasiparticle, with p = cjci —C2C2, f = 
{c[c2 + C2C1), where the c's are quasiparticle creation and destruction operators. The factor- 
ization or the semiclassical approximation (SCA) consists of assuming, equivalently, that the 
oscillator operators behave classically or that products of quasiparticle-oscillator operators 
can be factorized. 

We compare the SCA with the fully quantum mechanical results first by computing the 
polaron binding energies. This is done easily in the strong-coupling limit by freezing the 
quasiparticle hopping dynamics. We note first that, in the harmonic oscillator limit, the 
polaron binding energy is proportional to g^ whereas, in the opposite limit of the infinite 
square well, the width of the well remains finite and the interaction produces a lowering 
of energy that is proportional to g. This cross-over behaviour becomes evident from the 
full quantum- mechanical calculations as shown in Fig. 1. Plotted in the main figure is the 
binding energy (in arbitrary units) as a function of g. In the inset, the same quantities 
are plotted on a logarithmic scale. The bold lines indicate the fully quantum-mechanical 
calculation and the dashed lines indicate the results of the SCA. In all the cases, the oscillator 
frequency ujq has been kept fixed, and A varied allowing the oscillator to pass smoothly from 
the box (A = 1) limit to the harmonic (A -^ 00) limit. Two key results are evident in Fig. 1. 
First, the SCA results agree with the exact ones only in the harmonic oscillator limit; in the 
square-well limit, the departure becomes quite drastic. Second, (see inset), when the system 
is in the box-limit (A = 1), the fully quantum system behaves harmonically for small g but 
exhibits a cross-over for larger values of g showing the true box-limit slope of unity. 

We also calculate the overlap between the adiabatically displaced ground-state wavefunc- 
tions. This overlap, basically the Huang-Rhys factor, governs the polaronic tunneling rate 
in the strong-coupling limit. One knows, for instance, that in the harmonic oscillator limit, 
the tunneling rate is proportional to e^^ . In Fig. 2, we plot the overlap factor (logarithmi- 
cally) as a function of g for the box limit, i.e., A = 1 (dashed line), and for the harmonic 
oscillator limit, i.e., A ^ 00 (solid line). The quadratic dependence is clearly seen for the 



latter. However, for A = 1, the rise is sublinear, showing that the dependence of the overlap 
factor on g is much weaker. 

III. HEISENBERG EQUATIONS 

We discuss in this section the temporal evolution of the system, and show how the SCA 
differs from the fully quantum mechanical treatment. The equations of motion corresponding 
to the Hamiltonian @) can be written as 

p = 2Vq (6a) 

^ = -2Vp + 2uJog{X + l/2y/^fsmz (6b) 

f = -2uJog{X + 1/2)^/2^ sin z (6c) 

T^z (6d) 



(A + 1/2) 



Tiz = — u;o TT —7777- tan z sec^ z + g{\ + 1/2)^' ^p cos z \ , (6e) 

where q = i(c^C2 — clci) and the quasiparticle operators p, g, f cyclically satisfy the commu- 
tation relations, [p,q] = 2if, and [i,vf2] = i. As stated earlier, the SCA consists in assuming 
the oscillator operators z, ifz to be c-numbers. In the temporal analysis, we compare the re- 
sults of such an approximation with those given by the full quantum evolution described by 
Eqs. (H). We plot in Figs. 3-5, the evolution of the population difference of the quasiparticle 
between the two sites p{t) as a function of dimensionless time Vt. In all our calculations, 
the initial condition used for the quantum system is the ground state of the quasiparticle- 
oscillator system projected onto the one-site localized part of the Hilbert space, such that 
(p)(0) = 1. The initial condition used for the SCA calculation is p(0) = 1, i(0) = 'kziS^) = 0. 
In all the plots, the solid line indicates the full quantum evolution and the dashed line indi- 
cates the evolution due to the SCA. The polaron binding energy has been kept constant to 
facilitate comparison. This value (which we take to be 1.5 V) equals g'^uJo/2 in the harmonic 
oscillator limit. In Fig. 3, A = 200, the oscillator potential is essentially that of the harmonic 
oscillator potential. The oscillator energy ujq takes on the values 10 V, V, Q.IV in (a),(b), and 

6 



(c) respectively. As discussed elsewherea, whereas the SCA shows self-trapping for all the 
values of the oscillator frequency, the full quantum evolution differs substantially, except in 
the limit of low oscillator energy ujq = 0.11^. In this limit, for short times, the full quantum 
evolution and SCA agree in that both show self-trapping, with nearly the same average 
value of self-trapping and oscillation frequency. However, the quantum evolution shows a 
considerably richer structure involving collapses and revivals. At much longer times, the 
dressed quasiparticle tunnels from one site to the other. This is evident in Fig. 3 (b). When 
A = 10, (Fig. 4), the potential is more ' square- well'-like and some departures, especially in 
the small oscillator frequency regime (Fig. 4 (b)) are visible. For instance, the 'silent runs' 
separating the collapse and revival sequence, are less quiescent, and the agreement between 
the SCA and the full quantum evolution is slightly worse. In Fig. 5, we take A = 1, wherein 
the potential is essentially the infinite square-well. Whereas the agreement between the SCA 
and the full quantum evolution is best when the oscillator energy is least, ujq = O.IV^, (Fig. 5 
(c)), the agreement is far worse than for the harmonic potential (Fig. 5 (c)). Further, the 
'silent runs' are barely noticeable, with the collapses and revivals intruding into each other. 
A key time-scale in the temporal evolution is the one associated with the polaronic 
tunneling between the two sites. Since this time scale is intimately connected with the 
polaron binding energy, we plot in Fig. 6, the logarithm of the tunneling time as a function 
of the binding energy. The solid line denotes the harmonic oscillator limit, A ^ oo, whereas 
the dashed line denotes the infinite square- well limit, A = 1. Note that for small values of 
the binding energy, the time-scales for both cases is only weakly dependent on the energy. 
However, for larger coupling (binding energy), both show a clear linear dependence (albeit 
with different slope). This clearly indicates that even for the box-like potential, the polaronic 
tunneling time scale is exponentially dependent on the binding energy. While well-known 
for harmonic polarons, this exponential dependence constitutes an important new result for 
anharmonic polarons emerging from the present analysis. 



IV. SUMMARY 

By analyzing the dynamics and energetics of a quasiparticle interacting with a tunably 
anharmonic oscillator, specifically described by a Poschl- Teller potential, we find that, in 
the limit of strong coupling between the quasiparticle and the oscillator, selftrapping is 
robust and persists for strong anharmonicities, with the polaron tunneling time scale being 
exponentially dependent on the polaron binding energy, a feature that has been earlier 
known to be true for harmonic polarons. We further find that the full quantum result 
agrees with the predictions of the semiclassical approximation only in this strong-coupling, 
low-frequency regime, in agreement with earlier findings for harmonic polarons. 

One of us (VMK) acknowledges the financial support of the National Science Foundation 
under grant no. DMR-9614848, and of the Los Alamos National Laboratory under grant 
no. 0409J0004-3P. 
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FIGURES 
FIG. 1. Polaron binding energy as a function of the quasiparticle-oscillator coupling constant 

g. The inset shows the same quantities on a logarithmic scale. The solid line indicates the quantum 

mechanical result whereas the dashed line indicates the result of the semiclassical approximation 

(SCA). 

FIG. 2. The overlap of the adiabatically displaced groundstate wavefunctions plotted logarith- 
mically as a function of g for box limit, A = 1 (dashed line) and harmonic oscillator limit, A ^ oo 
(solid line). 

FIG. 3. The evolution of the quasiparticle probability difference p{t) as a function of dimen- 
sionless time Vt. In all the figures, the polaron binding energy has been kept fixed at 1.5 V and 
A ^ oo. The solid line denotes the fully quantum mechanical result and the dashed line denotes 
the result of the SCA. In (a), loq = lOV, (b), loq = V, (c)= loq = O.IV. 

FIG. 4. Same as Fig. 3, with A = 10. 

FIG. 5. Same as Fig. 3, with A = 1. 

FIG. 6. The polaron tunneling time scale plotted logarithmically as a function of the polaron 
binding energy, both computed without making the SCA. The solid line indicates A ^ cx) (harmonic 
oscillator limit) and the dashed line indicates A = 1 (box limit). 
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